The role of string-like, supramolecular assemblies in reentrant supernematic liquid 

crystals 
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Using a combination of isothermal-isobaric Monte Carlo and microcanonical molecular dynamics 
we investigate the relation between structure and self-diffusion in various phases of a model liquid 
crystal using the Gay-Berne-Kihara potential. These molecules are confined to a mesoscopic slit- 
pore with atomically smooth substrate surfaces. As reported recently [see M. G. Mazza et at, 
Phys. Rev. Lett. 105, 227802 (2010)], a reentrant nematic (RN) phase may form at sufficiently 
high pressures/densities. This phase is characterized by a high degree of nematic order and a 
substantially enhanced self-diffusivity in the direction of the director n which exceeds that of the 
lower-density nematic and an intermittent smectic A phase by about an order of magnitude. Here 
we demonstrate that the unique transport behavior in the RN phase may be linked to a confinement- 
induced packing effect which causes the formation of supramolecular, string-like conformations. The 
strings consist of several individual molecules that are capable of travelling in the direction of n as 
individual "trains" consisting of chains of molecular "cars". Individual trains run in parallel and 
may pass each other at sufficiently high pressures. 

PACS numbers: 61.30.Hn,66.10.C-,64.70.M-,61.30.Gd 



I. INTRODUCTION 



A thermodynamic phase is called "reentrant" when an 
order parameter varies nonmonotonicaly with the ther- 
modynamic field driving the transition. Reentrancy is 
ubiquitous in the physics of thermal many-particle sys- 
tems. It arises under quite disparate physical condi- 
tions encountered, for example, in quantum gases [ij, 
two-dimensional charged colloids Q , or relativistic scalar 
field models 13. As far as soft matter is concerned reen- 
trancy has been reported for self-assembled supramolec- 
ular structures j3j @i wetting phenomena at oleophilic 
surfaces 0], and novel discotic and calamitic liquid crys- 
tals 0. In fact, since the first observation of reentrant 
nematic (RN) phases in a seminal paper by Cladis Q 
reentrancy in liquid crystals seems to have received most 
of the attention. Liquid crystals offer rich and com- 
plex phase diagrams due to the characteristic interplay 
of translational and rotational degrees of freedom. For 
example, reentrant phase transitions have been reported 
for the isotropic phase of mixtures of discotic liquid crys- 
tals 0, the ferroelectric transition in smectic C phases 
, the cholesteric-to-blue phase transition in chiral liq- 
uid crystals [ll|, and for nematic (N) phases [l2| . 

After the work of Cladis @ on binary mixtures, and 
later in pure compounds at high pressure and in the su- 
percooled state [l^l, Hardouin et al. [3] showed the oc- 
currence of RN phases also for pure compounds at atmo- 
spheric pressure. For comprehensive reviews the reader 
should consult Refs. [l^ and [l^ Whereas most earlier 
work on reentrancy of phase transitions in liquid crys- 
tals is experimental in nature comparatively little atten- 



tion has been paid to this fascinating phenomenon from 
a theorist's point of view. Early theoretical approaches 
have mainly considered dimer models and frustration ef- 
fects |15l - [l8| . Ferrarini et al. consider the role of as- 
sociation and isomerization equilibria for the formation 
of RN phases within a mean-field- type of theory [l9j . 
The most recent theoretical study employs isothermal- 
isobaric and canonical ensemble Monte Carlo (MC) sim- 
ulations to investigate the nature of the smectic A(smA)- 
RN phase transition for a bulk system of hard ellipsoids 
with square-well attraction [2^ where earlier theoreti- 
cal studies are briefly reviewed, too. Unfortunately, the 
model employed in Ref. is somewhat artificial in as- 
suming that the ellipsodal molecules are always oriented 
in a perfectly parallel fashion such that all rotational de- 
grees of freedom are "frozen" irrespective of the thermo- 
dynamic conditions. Therefore, this study seems only 
of limited use to elucidate properties of RN phases at a 
molecular level. Moreover, the authors do not consider 
dynamic features of RN phases. 

Despite the variety of systems and thermodynamic 
conditions under which reentrancy in liquid-crystalline 
materials arises comparatively little attention has been 
paid to the dynamics of reentrant phases. To the best 
of our knowledge the only study of the dynamics of 
reentrant liquid-crystalline phases has recently been con- 
ducted by ourselves and focuses on self-diffusion in the 
RN phase as a specific case [2l|. In that work we con- 
sidered rod-like mesogens where the interaction between 
a pair of rods takes into account both a realistic shape 
of the molecules and the orientation dependence of their 
interaction. Unlike dc Miguel and Martin del Rio [2^ 
we allowed our molecules to rotate as freely as the spe- 
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cific thermodynamic conditions permit. We could then 
demonstrate [2l| that on account of the high degree of 
nematic order diffusion of molecules in the direction of 
the nematic director may be enhanced by about an order 
of magnitude over that characteristic of lower-density N 
or smA phases. The combination of a high degree of ne- 
matic order with a substantially enhanced sclf-diffusivity 
in the direction of the nematic director prompted us to 
refer to liquid crystals in the RN phase as "supernemat- 
ics" [2l| . Enhanced self-diffusivity in the RN phase bears 
a striking similarity to "levitation" in zeolites. Experi- 
ments have demonstrated the self-diffusion constant to 
pass through a maximum when the size of a diffusant 
molecule is comparable to the pore size, that is enhanced 
diffusivity is observed as a result of severe confinement to 
tiny voids in the host material (see, for example, Ref . [22h . 

For liquid crystalline materials we could explain this 
unusually enhanced self-diffusion quantitatively in terms 
of a severely reduced rotational configurational en- 
tropy [U . The relation between self-diffusion and struc- 
tural features of the RN phase has not been explored. 
In fact, as we shall show in this work the RN phase is 
characterized by unique conformations of entire groups 
of its molecules to which we refer as "strings" . Strings 
may be thought of as train-like arrangements of several 
molecules that stay together as molecular "cars" and dif- 
fuse as a supramolecular entity. As we argue below it 
is the fairly high pressure that keeps the trains together 
and which allows them to overcome the attractive in- 
teractions between molecules pertaining to neighboring 
trains. From the experimental side there also exists some 
evidence suggesting unique dynamic features of the RN 
phase. For example, distinct differences in the molecular 
dynamics in the N and RN phases can be concluded from 
corresponding changes in the nuclear magnetic resonance 
(NMR) relaxation times reported in Rcfs. [23r.26. How- 
ever, a more detailed account of the present experimental 
state of the art is postponed until Sec. IIVI 

Last but not least we note that an improved under- 
standing of the dynamics of liquid crystals is not only 
of academic interest but may also be important from a 
purely practical point of view. For example, by pumping 
liquid crystalline material through vesiscles, that maybe 
perceived as confined geometries, a spider is capable to 
spin silk with unsual materials properties that result es- 
sentially from the transport of the liquid crystal. Mate- 
rials with properties comparable to spider silk cannot be 
produced in the laboratory to date (see Fig. 2 in Ref. [l^- 

The remainder of our manuscript is organized as fol- 
lows. Our model liquid crystal is introduced in Sec. |TT] 
whereas Sec. IIIII is given to a presentation of our find- 
ings. In Sec. IIVI we summarize our results, discuss their 
physical origin, and put them in a broader experimental- 
theoretical context. Specifically, we rationalize the re- 
lation between the formation of supramolecular strings 
and enhanced self-diffusion in the RN phase. Details 
about the theoretical background of our simulations can 
be found in Appendix [Xj 



II. SIMULATIONAL DETAILS 
A. Model potentials 




FIG. 1: Cartoon of a pair of spherocylinders of length L -\- a 
where Ui and Uj describe the orientation of molecules i and j 
of that pair in a space-fixed Cartesian coordinate system; Vij 
is their center-of-mass distance vector which can be decom- 
posed into components (anti-)parallel (t'I'j) and perpendicular 
(r^) with Ui, respectively, as indicated. 

As in [U we employ a fluid composed of N sphero- 
cylinders of length L and diameter a capped at both 
ends by a hemisphere of the same diameter 2R = a (see 
Fig. [T]) . For this system the Hamiltonian may be cast as 
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where P — {pi,p2, ■ ■ . ,Pn) are the linear momenta 
conjugate to the set of center-of-mass positions R = 
(ri,r2, . . . ,rAr) and T = (71,72, . . ■ ,7Af) are the Eulcr 
angles specifying the orientation of each molecule. For 
the special case of spherocylinders with two rotational 
degrees of freedom, 7^ = {9i,(pi) is a two-dimensional 
vector composed of two of the three Euler angles. In 
Eq. ([l} , Cia is the component of the angular momentum 
referred to the body- fixed principal axis a of molecule i, 
m is the molecular mass, and 
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is the moment of inertia of a sphero cylinder. 

In a computer simulation [regardless of whether Monte 
Carlo (MC) or molecular dynamics (MD) is employed] 
the key quantity is the configurational potential energy 
U{R,U) of the system [see, for example, Eqs. (|A2p . 
(|A5b[) . or (|A7|) ] where U = (mi, M2, • ■ • ,Un) is a set of 
unit vectors specifying the orientation of the molecules 
in a three-dimensional space-fixed Cartesian coordinate 
system (see Fig.[T]). Elements of the set U are related to 
corresponding elements of T via 



dui = sin 6id6i 



1, 
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In the general case of a confined fluid U{R, U) can be 
decomposed into a fluid-fluid (ff) and a fluid-substrate 
(fs) contribution. Assuming pairwise additive interac- 
tions the former can be cast as 
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whereas the latter may be expressed as 
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k=l i=l 
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In Eqs. (jl]) and (O, and itfs represent the specific 
model potential adopted to describe the intermolccular 
interaction between two rodlike molecules and between 
a rodlike molecule and a planar, structureless solid sub- 
strate, respectively. In this work we adopt the so-called 
Gay-Berne-Kihara (GBK) model In the GBK model 
the interaction between a pair of spherocylinders depends 
on the relative molecular orientation through the func- 
tion d^{rij ,Ui,Uj) which depends on the centcr-of-mass 
distance vector = rt ~ rj and the orientations of 
molecules i and j such that c?°J is actually the minimum 
distance between that pair of molecules. More specifi- 
cally. 
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where r = r/r and r = |r|. In Eq. ^ the function 



esirij,Ui,Uj) = gff <^ I - — - 
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describes the orientation dependence of the anisotropy 
strength where the parameters x ^-nd x' given by 
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In these latter two expressions parameters k, = L + 1 {L 
in units of a) and k' may be thought of as the aspect ra- 
tio of a spherocylinder and the interaction strength for a 
side-side relative to an end-end configuration of a pair of 
spherocylinders, respectively. This can be seen from vari- 
ous plots in Fig.[2]indicating that the potential minimum 
depends strongly on the molecular orientation such that 
a side-side configuration of a molecular pair is the ener- 
getically most favored one. One also notices from Fig. [5] 




FIG. 2: (Color online) Plots of tiff as functions of minimum 
distance dfj{rij,Ui,Uj) between a pair of spherocylinders in 



units of eg and a, respectively, and for various relative orien- 
tations such as side-side (Fij • Ui = rij ■ Uj = 0, Ui ■ Uj = 1) 
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T-shaped {rij ■ Ui = 0, fi. 
crossed {vij ■ Ui = G, Vij 
Fig.©. 
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that both the location of the potential minimum as well 
as the range of distances over which becomes repul- 
sive depend very little on the intermolecular orientation 
as one would expect on account of the spherocylindrical 
shape of the molecules. 

In a similar spirit, d^{zik, Ui) introduced in Eq. ([5]) is 
the minimum distance between a spherocylinder with ori- 
entation Ui whose center of mass is at a vertical distance 
Zik = Zi± Sz/2 from the lower (k = 1) and upper (fc = 2) 
substrate located at — Sz/2 and -f-Sz/2 along the z-axis of 
a space-fixed Cartesian coordinate system, respectively. 
Assuming the substrates to be planar, d"^ can easily be 
determined as the smaller of the two distances of both 
end points of a spherocylinder from the substrate plane. 
The corresponding function d™(ry, tii, Mj) is more com- 
plex but can be computed numerically u sing an efficient 
algorithm proposed by Vega and Lago [23|- Following 
l2ll we model the fluid-substrate interaction via 
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(9) 



where the parameter efs controls the strength of interac- 
tion similar to eff in Eq. ([7]) and ps ~ 2/£^ is the areal 
density of the substrate where (/a = is the lattice 

constant of a single layer of atoms arranged according to 
the (100) plane of the face-centered cubic lattice. The 
diameter a of these substrate atoms is taken to be the 
same as the diameter of a spherocylinder of the confined 
fluid phase. To obtain Eq. © we assume the interac- 
tion between the point on the spherocylinder closest to 
the nearest substrate atom to be described by a conven- 
tional Lennard- Jones potential. Holding both the center- 
of-mass position and the orientation Ui fixed one then 
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averages the interaction with the substrate atom over the 
entire substrate area. Notice that this procedure gives 
rise to a prefactor | in front of the repulsive term in 
Eq. © which we have dehberately omitted for conve- 
nience. 

In Eq. (HJ, < g{u) < 1 is the so-caUed "anchor- 
ing function" . It permits to discriminate energetically 
different orientations of a molecule with respect to the 
substrate plane. In this paper we employ 

g{u) = {u ■ e^f + {u ■ eyf (10) 

where Cq is a unit vector pointing along the a-axis of 
a space-fixed Cartesian coordinate system. Hence, any 
molecular arrangement parallel with the substrate plane 
is energetically favored whereas a homcotropic alignment 
of a molecule (u || e.^) receives an energy penalty by 
"switching off" the fluid-substrate attraction altogether. 
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FIG. 3: (Color online) Plots of the equipotential contour plot 
'ifi/efi for the GBK (bottom), K (top left), and GB (top 
right) models (see text); color code shown alongside. The 
area shaded in white is defined through the inequaility itft > 
therefore approximately representing the molecular shape. In 
generating the plots we assumed a pair of molecules to be 
located in the x-y plane with a perfectly parallel orientation 
(i.e., Ui-Uj = 1) as a special case (see also Fig.[l|. Both x and 
y are given in units of a. Note that the GBK model combines 
the molecular shape of the K model with the anisotropy of 
the GB model. 

At this stage a couple of additional comments apply. 
For example, the reader may notice that Eqs. ([7]) and @ 
are identical with those describing the orientation depen- 
dence of the interaction strength in the well-known Gay- 
Berne (GB) model of liquid crystals [13]. However, unlike 
in the GBK model a liquid crystal molecule is perceived 



as an ellipsoid of revolution in the GB model. This is ef- 
fected by replacing the function dfj{rij,Ui, Uj) in Eq. ^ 
by one that accounts properly for the shape of a molecule 
in the GB model [see, for example, Eq. (4) of Ref. [sij . 
The GBK model is also closely related to the so-called 
Kihara (K) model ^3 of a liquid crystal [s^. In this 
latter model a molecule is again perceived as a sphero- 
cylinder but the interaction strength between a pair of 
molecules depends only on the minimum distance but 
not on the relative orientation of the molecules of a pair. 
This is achieved by setting k' = 1 in Eq. (jSbp such that 
eff{rij,Ui,Uj) = eg in Eq. ([7]). The shapes of molecules 
in the K, GB, and GBK models resulting from these ma- 
nipulations are illustrated by the contour plots in Fig. [3l 



B. Numerical details 

The simulations (both MC and MD) to be presented 
below are based upon systems containing N = 1500 
molecules. Quantities of interest will be expressed in 
customary dimensionless (i.e., "reduced") units. For ex- 
ample, length will be expressed in units of cr, energy in 
units of eg, temperature in units of eg/kB, time in units 
of (cr^m/eff )^/^ using m ~ 1, and pressure P|| in units of 
k-QT/a^ where P\\ = ^(Pxx + ^yy) is related to diagonal 
components of the pressure tensor P acting in the x-y 
plane. 

We employ MC simulations in an isothcrmal-isobaric 
ensemble where we use iV, T, P|| , and s^q input param- 
eters. In all the simulations we fix Szo = 19 and consider 
T ~ 4.0 and 6.0. To reduce the computational cost in- 
teractions between fluid molecules are cut off beyond a 
minimum distance = 3. In addition, we employ a 
conventional Verlct neighborlist where molecules up to a 
minimum distance of 8 are included as neighbors [3J] . In 
the X- and y-directions periodic boundary conditions are 
applied. 

In MC we employ the algorithm of Schoen [s^ with 
one important modification. Rather than changing the 
area A of the simulation box in the sense of a similarity 
transformation we allow the individual side lengths 
{a = X, y) to vary independently such that the shape 
of the simulation cell may change during the course of 
a simulation. This is important to preserve the in-plane 
isotropy of P especially in highly ordered N, smA, or 
RN phases. Our runs start from a sufficiently low pres- 
sure in the isotropic phase. In the initial configuration 
molecules are perfectly aligned and then allowed to relax 
to the isotropic state. The last configuration of this run 
is then used as a starting configuration of the next run 
at a slightly higher P||. We continue with such a com- 
pression sequence until the last state point of the highest 
pressure of interest has been reached. We refer to a MC 
step as N attempted displacements or rotations and one 
attempted change of Sq,. Displacements and rotations 
as well as changes for a = x or a = y are performed 
with equal probability. Our runs are based upon 3x10^ 
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equilibration steps followed by 2 x 10^ production steps 
for thermodynamic states sufficiently far away from any 
phase transition; around these transitions the number of 
production steps was enlarged to 7 x 10^. However, in 
changing A as described care has to be taken that none of 
the instantaneous side lengths Sa becomes shorter than 
twice the neighbor-list cutoff given above. 

MC simulations arc mainly employed to generate suit- 
able starting configurations for the subsequent MD sim- 
ulations. In addition, MC is used to independently verify 
the correctness of the MD simulations through a compar- 
ison of equilibrium properties accessible to both types of 
simulations. In MD we employ the velocity Verlet al- 
gorithm for linear molecules in the implementation sug- 
gested by Ilnytskyi and Wilson [s^. The iterative so- 
lution of the equations of motion is initiated by using 
the last configuration of the MC simulation at the same 
state point as a starting configuration where velocities are 
taken at random from a Maxwell-Boltzmann distribution 
at the desired temperature. This starting configuration is 
equilibrated for another 10^ time steps of St = 0.02 using 
a simple velocity rescaling thermostat. The equilibrated 
system is then monitored for another 2 x 10'' time steps 
with no applied thermostat. During this part of the run 
the integration time step is reduced to St = 10^'^. 



III. RESULTS 

A. Structure of ordered liquid-crystalline phases 

To characterize the structure of N, smA, and RN 
phases we begin by introducing suitably defined order 
parameters. For nematic phases this is accomplished via 
the so-called alignment tensor defined as [s^l 



Q 



1 " 



2N 
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(11) 



where "(g)" denotes the direct (i.e., dyadic) product and 
1 is the unit tensor. Hence, Q is a real, symmetric, and 
traceless second-rank tensor which can be represented by 
a 3 X 3 matrix. For a given configuration (i?, U), Q can 
be diagonalized in the basis of its three eigenvectors n_, 
no, and ri_|_ with associated eigenvalues A_, Aq, and A+. 
Conventionally jssj . one then defines the Maier-Saupe 
[1^,1131 nematic order parameter S = (A+) where angular 
brackets denote an ensemble (MC) or time average (MD). 
On account of its definition, 5 = in the isotropic phase 
if the orientation of the molecules is perfectly random 
whereas S* = 1 if the molecules arc perfectly aligned with 
the director n = in the nematic phase. 

Even though the nematic phase is characterized by 
a substantial degree of orientational order the distribu- 
tion of center-of-mass positions exhibits the typical short- 
range order characteristic of ordinary fiuids where the 
typical correlation lengths arc comparable to the range 
of intcrmolecular interactions. This changes in the smA 




FIG. 4: (Color online) Plots of nematic S (open symbols) 
and smectic order parameter A (filled symbols) as functions 
of pressure Py for various liquid-crystalline phases, T — 4.0 
(circles), and T — 6.0 (squares). I (blue symbols), N (red sym- 
bols), smA (green symbols), and RN phase (black symbols). 
The same color code for various liquid-crystalline phases is 
used throughout this work. 



phase where molecules do not only align with n but de- 
velop long-range positional order, that is they form indi- 
vidual layers in the direction of n. However, positional 
order remains short-range in directions orthogonal to n 
such that within a given layer the structure remains fluid- 
like. Because of the layering of the fluid in the direction 
of n a suitable quantitative measure of smectic order is 
given by the leading coefficient of the Fourier transform 
of the local density p{r ■ n) 



A = 



1 

N 



N 



^exp 



27Ti {Vi 



n 



(12) 



where d is the spacing between adjacent smA layers. If 
these layers were ideal, d = k (see Sec. Ill Ap . However, 
in practice d > k on account of thermal fiuctuations. 
Hence, in each configuration {R,U), d is adjusted such 
as to maximize A. From Eq. (jl2p it is also apparent that 
As [0, 1]. In the nematic phase the superposition of the 
complex exponential functions in Eq. (fT2|) causes A « 
because in the direction of n a periodic structural feature 
of peridodicity d does not exist; by a similar token Awl 
in the smA phase. 

Plots of both S and A arc presented in Fig. 2] for tem- 
peratures T = 4.0 and 6.0. As one can see from the 
plots, both S and A are small at low pressures indicat- 
ing the presence of the I phase. Ideally, 5* = in the 
I phase but in practice attains small positive values of 
about 0.2. This is partly due to the substrates which 
induce a certain degree of nematic order in those por- 
tions of the confined liquid crystal that are closest to 
either substrate. However, a nonvanishing value of S in 
the isotropic phase is also due to a small but significant 
finite-size effect that is well-known and has been studied 



quantitatively in the literature [s^ HTj . Here we adopt 
a value S" ~ 0.4 as a sensible heuristic definition of the 
threshold above which the N phase exists in agreement 
with the mean- field theory of Maier and Saupe [s^ . 
From Fig. |4]it therefore appears that the N phase forms 
somewhere above P|| ~ 1.0 (T = 4.0) and 1.2 (T = 6.0), 
respectively. 

Focusing on the lower temperature first, S increases in 
the N phase with pressure up to P|| ~ 1.7 indicating the 
increase of nematic order. Over the same pressure range 
there are no smectic layers because A < 0.2 remains rel- 
atively small. However, at P|| ~ 1.8, smectic layers arc 
forming as reflected by A ~ 0.8. Simultaneously, S keeps 
increasing as well because in the smA phase molecules 
align themselves even better with n. The smA phase re- 
mains stable until P|| ~ 2.2 has been attained. Beyond 
this pressure S remains nearly constant at a high value of 
about 1.0 whereas A drops rapidly until a residual small 
value of about A ~ 0.1 is assumed which remains con- 
stant at all higher pressures considered. The drop in A 
with P|| clearly reflects the disappearance of smectic lay- 
ers present at lower P||. Because 5' ~ 1.0 over the same 
pressure range we conclude that at sufficiently high pres- 
sures the structure of the confined fluid becomes nematic 
again in a reentrant fashion. To distinguish the high- 
pressure from the low-pressure nematic we refer to the 
former as RN because N and RN phases are seperated 
by an intermittent smA phase. 

At T = 6.0 these general features of the order- 
parameter plots in Fig.Ulprevail except that A never rises 
above the residual value of about 0.1. Hence, the smA 
phase does not form at this temperature. Nevertheless, 
we notice a small, step-like increase at P|| ~ 2.2 in the 
plot of A. At this and all larger pressures considered the 
nematic order is high and increases even further reflected 
by a monotonic increase of 5* toward its limiting value 
1.0. For P|| > 2.2 the confined fluid exhibits structural 
features of the RN phase (see below). We shall therefore 
keep the acronym despite the missing intermittent smA 
phase. 

To get a better feel for the structure of the confined 
fiuid in the N, smA, and RN phases we present con- 
figuration "snapshots" from our simulations in Fig. m 
The plot in Fig. [SJa) illustrates the structure of the N 
phase: one clearly recognizes a preferred net orientation 
of the molecules but no significant positional order of the 
center-of-mass distribution. In Fig. [SJb) the layer struc- 
ture characteristic of the smA phase is visible. These 
layers disappear in the RN phase as shown in Fig. [SJc). 
However, a large degree of orientational order is preserved 
as one realizes from this latter plot. 

A more quantitative analysis of the structure of con- 
fined liquid crystals in the various ordered phases is en- 
abled by plots of the local density 

= (13, 
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FIG. 5: (Color online) "Snapshots" of characteristic configu- 
rations at r = 4.0; (a) Py = 1.6 (N), (b) = 2.1 (smA ), (c) 
P|| = 2.9 (RN) (see Fig. |4|. Molecules in green are aligned 
with fi whereas rectangular areas in dark green represent the 
solid substrates. 



where S {z ~ zi) denotes the Dirac iS-function, TV (z) is 
the number of molecules with their centers of mass lo- 
cated within a interval of width (5z = 0.1, and A = SxSy 
is the instantaneous area of the simulation cell of side 
length Sa [a = x, y) in the a-direction. Hence, p (z) is a 
measure of the probability of finding the center of mass 
at a specific position z. Plots in Fig. [5] illustrate the 
local structure of the confined liquid crystal in the direc- 
tion of the substrate normal (i.e., perpendicular to the 
director n). As is well-known from simple fiuids (42j . in 
confinement p(z) is an oscillatory function of position suf- 
ficiently close to the substrate surfaces. The oscillations 
indicate that molecules arrange their centers of mass in 
individual layers parallel with the substrate. The loca- 
tion of a molecular layer is therefore indicated by a peak 
in p (z). In the immediate vicinity of the solid substrate 
p (z) because of the diverging fiuid-substrate repul- 
sion as \z\ — >■ Sz/2. 

These general features may be seen from the plots in 
Fig. [B] One notices that in the N phase oscillations 
in p [z) are restriczed to the immediate vicinity of the 
solid substrates [see Fig. HJa)]. The oscillations decay 
exponentially as \z\ — >■ with a correlation length de- 
termined by bulk properties [43| whereas the spacing 
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FIG. 6: (Color online) Local density p (z) as a function of z 
for T = 4.0; (a) Py = 1.7, (b) = 2.2, (c) P,, = 2.3 (see 
Fig. 13] for color code). 



between neighboring layers reflects the molecular shape 
j4^ . In the N phase a bulk-like region of considerable 
width exists in which p (z) is constant as the plot in 
Fig-inja) reveals. On the contrary, layering is much more 
pronounced in the smA phase as a comparison of plots 
in Figs. [6l^a) and [6l[b) indicates: Peaks in Fig. El^b) are 
taller and because a bulk-like region is absent the en- 
tire confined fluid is organized in layers. However, the 
layered structure becomes weaker as one approaches the 
center of the slit pore as reflected by the decay of p (z) in 
Fig-El^b). Finally, in the RN phase layering is even more 
pronounced as the plot in Fig. |6l^c) shows. Moreover, 
plots in Fig. El^b) and Hfc) reveal another subtle struc- 
tural change. Counting maxima in the plots of p {z) one 
realizes that one additional layer of molecules is present 
in the RN compared with the smA phase. The presence 
of this additional layer points to a more efficient packing 
of molecules as far as the arrangement of their centers 
of mass is concerned. A more efficient packing becomes 
possible because the layered structure characteristic of 
the smA phase [see Fig. E^b)] is lost in the RN phase [see 
Fig.EIc)]. 

Next, at T = 6 one realizes from Fig. [7] a similar 
formation of an additional molecular layer but without 
occurrence of an intermittent smectic phase. Plots in 
Figs. [7l[a) show that for two state points in the nematic 
phase an increase in pressure only causes the layers al- 
ready present to become more distinct. For example, 
p{z) drops to zero in regions between individual layers 



FIG. 7: As Fig. El but for T = 6.0; (a) Py = 2.3 ( ) and 2.4 

(---), (b) P|| = 2.5, (c) P|| = 2.6 (see Fig. [4] for color code). 
Note that all parts of the figure refer to RN states. Inset is 
a magnification of curves for P|| = 2.4, 2.5, and 2.6 around 

2 = 0. 



at the higher pressure whereas p (z) remains nonzero over 
the corresponding ranges at the lower pressure. The plot 
in Fig. [7|b) corresponds to a state point where the con- 
fined fluid undergoes the subtle reorganization that will 
eventually lead to the growth of an additional layer. As 
one approaches the midplanc located at z = individ- 
ual peaks in p (z) first broaden and eventually split up 
into two unresolved peaks. At an even higher pressure 
considered in Fig. we realize that the fluid has com- 
pleted its structural reorganization and a new fluid layer 
emerges as we can see by comparing this plot with its 
counterparts shown in Fig. [TJ^a). Moreover, visual in- 
spection of a characteristic configuration [see Fig. [SJc)] 
reveals that in the RN phase molecules appear in string- 
like arrangements. 

In our terminology a "string" is a sequence of molecules 
arranged in space like individual cars of a train with their 
long axes all pointing in the same direction. The "snap- 
shot" of a typical configuration presented in Fig. [SJc) 
already indicates that these supramolccular strings are 
omnipresent in the RN phase. In mathematical terms an 
ideal string can therefore be defined through the relations 



Ui ■ U i 



= 1 

= Mo • r, 



1 



(14a) 
(14b) 



where, of course, one of the two scalar products in 
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a) 




FIG. 8: (Color online) Sketch of a pair of spherocylinders 
(see Fig. [T|. (a) Shifted side-side conformation, (b) As (a), 
but for a string-like conformation. The wedge-shaped areas 
corresponds to the solid angle Ui ■ Uj > S^ob b used in Eq. p6p 
which we exaggerate for illustrative purposes. 



Eq. (|14bp is redundant. In practice, however, small de- 
viations from perfect alignment of molecules in a string 
arise on account of thermal fluctuations. A good measure 
of the extent to which a pair of molecules is arranged in 
a string-like conformation approximately may be defined 
by decomposing the distance vector r^j between the cen- 
ters of mass of nearest-neighbor molecules i and j into 
components parallel and perpendicular to Ui through the 
expressions (sec Fig. [T]) 



u, ■ ri 



(15a) 
(15b) 



For a string-like configuration of a pair of molecules 



one anticipates rf. 



0. However, a second 



characteristic conformation characterized by relatively 
small values of r,4; is conceivable. In this conforma- 
tion molecules are arranged in a shifted side-side arrange- 
ment [see Fig. Hl^a)] . Plots of configurational snapshots 
in Fig. [5] clearly show that also unshifted side-side ar- 
rangements exist both in the N and in the RN phase. 
This is because they are energetically favored according 
to Fig. [3]Ja). However, they are also present to some ex- 
tent in other phases as snapshots in Figs. [5ja) and [5][c) 
indicate. To eliminate these unwanted unshifted side-side 
conformations in favor of only shifted ones it is prudent 
from an operational point of view to consider only the 
subset {r^} C {r,^} defined through 



< r,^ = K + <5« A r|i . > - + 1 



AUi ■ Uj > Scose} 



(16) 



where we arbitrarily take S^. = 0.2 and 6cose ~ 0.98. 
To see that this rather complex definition selects only 
those conformations sketched in Figs, ^a.) and EUb) one 
needs to realize that the first of the three conditions in 



J. 




FIG. 9: (Color online) Probability distribution V of (see 
Figs. [U and m Eqs. ^-^]. (a) T = 4.0 and P,, = 2.3 
(V), P|| = 2.6 (□), and Py = 2.8 (O); inset: Py = 0.9 (+), 
P|| = 1.4 (x), and Py = 1.9 (A), (b) T = 6.0 and Py = 2.2 
(o), Py = 2.4 (V), P|| = 2.5 (□), and Py = 2.9 (O); inset: 
Py — 1.1 {+) and P|| = 1.7 (x). For line and symbol color 
code see Fig. [l] 



Eq. (fT6|) selects nearest neighbors {j} of a given molecule 
i, the second condition favors head-tail-like conforma- 
tions, whereas the third one selects molecules with nearly 
perfectly aligned long molecular axes. It is then conve- 
nient to introduce the probability 



\l=l I 



(17) 



for a quantitative discussion of the formation of string- 
like structures where Vo is a normalization constant 
which we determine via 



dr^p(r-L) = 1 



(18) 



In the actual simulations we obtain 'P{r-^) as a histogram 
using 5r^ = 0.02 as the bin width. 
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Plots in Fig. |9]show V{r^) for different temperatures 
and pressures in the I, N, smA, and RN phase. At suf- 
ficiently low pressures it is evident from the insets in 
Figs. [5Ja) and ini[b) that V{r^) increases monotonically 
with and reaches a broad maximum located approxi- 
mately at ~ 1.25 irrespective of whether the I, N, or 
smA phase is considered. That V{r^) rises with r-^ may 
be ascribed to molecules in the I, N, and smA phase that 
arrange their long axes locally to a certain degree. This 
local order is solely caused by the intermolecular interac- 
tion potential which favors local side-side conformations 
of molecules [see Fig.[3][a)]. 

In the RN phase the form of 7'(r^) is completely dif- 
ferent. Here the plots in Fig. |9l[a) and El^b) exhibit two 
maxima. The largest one slightly above r-^ ~ 1.0 cor- 
responds to conformations in which neighboring parti- 
cles are arranged in the shifted side-side conformation 
depicted schematically in Fig. [5l^a). Note that these val- 
ues of T"*- are slightly below the position of the minimum 
of the interaction potential (see Fig. [2]). The secondary 
maximum of V{r-^) in the RN phase at < 0.5 cor- 
responds to string-like conformations schematically illus- 
trated by the plot in Fig.|5l^b). Were the string-like con- 
formation perfect one would anticipate a peak of 'P{r^) 
at = 0. However, on account of thermal fluctuations 
these ideal string-like conformations are only infrequently 
observed. For example, closer scrutiny would reveal that 
V{r^) 7^ at r-'- = indicating that a very small frac- 
tion of molecules does, in fact, assume an ideal string- 
like conformation. As P|| increases one notices that both 
peaks of V{r^) shift to progressively smaller . This 
fact reflects an arrangement of molecules in more perfect 
string-like conformations with increasing pressure as one 
would expect intuitively. 

One also notices from plots in Figs.lHl^a) and^l^b) that 
the curves seem to intersect in isolated points at char- 
acteristic values of r-^. These intersections may be in- 
terpreted as "isosbestic" points (derived from the Greek 
words isos: equal, the same and shestos: extinguishable) 
psj . a concept widely known in spectroscopy. The oc- 
currence of isosbestic points is relatively simple to grasp. 
We assume that V{r^) can be decomposed into a sum of 
three contributions: 

1. a probability distribution centered at r"*" « 1.0 
corresponding to a population of shifted side-side 
molecular conformations [see Fig. [5Ka)], 

2. a probability distribution centered at 7'^ ss 0.25 
corresponding to a population of shifted string-like 
molecular conformations [sec Fig. [8Kb)], 

3. and a probability distribution representing all other 
molecular conformations satisfying Eq. ([TS]) and 
embracing the entire range of r-^ . 

For sufflciently large P|| the distributions of [T] and [2] 
can be assumed to be well separated whereas |3] may be 
thought of as a broad background distribution. Inter- 
sections between distributions [T] and [3] or between [2] and 




FIG. 10: (Color online) (a) Persistence probability a (r) as a 
function of time r and T = 4.0; P\\ = 0.9 (o), P\\ = 1.7 (□), 
P|| = 2.1 (O), P|| = 2.3) (A), P|i = 2.8 {+). (b) as (a), but for 
T = 6.0 and Py = 1.1 (o), Py = 1.9 (□), Py = 2.4 (a), and 
P|| = 2.9 (+) (see Fig.glfor color code). 

[HI respectively, will give rise to isosbestic points. As P|| 
increases distributions [1] and [D become sharper and in- 
crease in height at the expense of the background dis- 
tribution [3] furthermore, maxima of both distributions [1] 
and [21 shift to lower values of indicating an increase 
in order. Importantly, the emergence of two wcll-deflned 
peaks in P(r^) is a fingerprint of the RN phase and could 
therefore be taken as computational evidence to detect 
RN phases even though its experimental significance may 
be limited. Incidentally, at low Py (i.e., in the I, N, or 
smA phase) the only important distributions are [T] and 
[21 because string-like conformations are statistically irrel- 
evant. Consequently, only one isosbestic point exists as 
shown in the insets of Fig. [5] 



B. Self-diffusivity in liquid crystalline phases 

The formation of string-like conformations has impor- 
tant ramifications as far as dynamic properties of the RN 
phase are concerned. To illustrate this we begin by con- 



10 



sidering a persistence probability 



air) 



J2 + ^) 



(19) 



defined as the correlation function of the angle between 
a pair of molecules which at the initial time t sat- 
isfy the conditions for a string-like conformation intro- 
duced in Eq. p^ . Subscript "i" is attached to the angu- 
lar brackets to emphasize that averaging involves a suf- 
ficiently large number of statistically uncorrelated time 
origins which implicitly invokes the principle of station- 
arity of temporal correlation functions in equilibrium sys- 
tems [4^ . In Eq. ([T9| "qij = Ui ■ Pij is the cosine of the 
angle between the instantaneous orientation of a refer- 
ence molecule i and the unit vector connecting its center 
of mass to that of another molecule j as time passes [see 
Fig- [I]- Because of this definition, rjij {t -I- r) is a measure 
of the relative orientation of a molecular pair at time t+T 
given that its orientation at time t is rjij (t). Therefore, 
cr (r) — > 1 as r — > in agreement with plots in Figs.fTOla) 
and llOf b) for various temperatures and pressures. At 
sufficiently low P|| , a (r) eventually decays with increas- 
ing time. The decay is delayed the more as the nature 
of the liquid-crystalline phase changes in the direction 
I^N^'SmA, that is with increasing P|| [see Fig. [TUT a)]. 

One also notices from Fig.fTUTa) that the decorrelation 
of relative orientations becomes very slow as soon as one 
enters the RN phase. This is already evident from the 
plot for P|| = 2.3 immediately after the smA-RN phase 
transition has occurred (see Fig. 2]). However, the cor- 
responding plot in Fig. [TITF a) shows that (t(t) eventually 
decays at long times. This is apparently not so at the 
higher pressure P|| = 2.8 in the RN phase. The rele- 
vant plot in Fig. [TUT a) shows only a weak decay of a (r) 
over an intermediate time range 1 ^ t < 10 and then 
remains nearly constant for longer times r > 10^. The 
shape of a (r) seems to refiect a two-stage decorrelation 
process involving those conformations represented by the 
bimodal distribution V{r-^) plotted in Fig.[9la). In fact, 
intuitively it seems sensible to assume that decorrelation 
during the intermediate time interval is associated with 
the orientational dynamics of shifted side-side conforma- 
tions depicted in Fig. [DJa). These conformations are ex- 
pected to be less stable than the string-like ones depicted 
in Fig. El^b). This notion is further corroborated by the 
slowing down of the dynamics at longer times with in- 
creasing P|| which causes the population of string-like 
conformations to increase according to plots in Figs. 1^1 
This qualitative picture remains unaltered at the higher 
T = 6.0 as parallel plots in Fig. [TUTb) show. However, at 
this temperature the smA phase is absent in agreement 
with our earlier analysis of various structural quantities. 

From the discussion of orientational relaxation one 
may get the impression that the RN phase is rather 
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FIG. 11: (Color onlme) (a) MSDs (Arjj (r)) as functions of 
time r and T — 4.0 and various liquid-crystalline phases at 
P|| = 0.9 (I) (o), P|| = 1.7 (N) (□), P|| = 2.1 (smA) (O), and 
P|l = 2.6 (RN) (A), (b) as (a), but for T = 6.0 and Py = 1.1 
(I) (o), P|| = 1.9 (N) (□), and Py = 2.4 (RN) (A) (see Fig. H 
for color code). 



sluggish with respect to its dynamics. While this is cer- 
tainly true for the rotational dynamics it turns out that 
mass transport is actually quite fast in the RN phase. A 
suitable measure of mass transport is the mean square 
displacement (MSD) of molecules. It measures the on- 
average net distance a molecule's center of mass is able 
to travel in a given amount of time. Because of the large 
degree of ncmatic order characteristic of the RN phase it 
seems sensible to introduce a specialized MSD by consid- 
ering displacements of molecules in the direction of their 
long axes. Specifically, we define r| = Ui ■ ri such that 
the associated MSD for the displacement of molecules in 
the direction of their long axes may be cast as 



Ar (r) ^ 



+ r)~r^ it) 



(20) 



Plots in Fig. [TT] illustrate the variation of the MSD 
with time on a double-logarithmic scale. At short times 
T < 10~^ MSDs for various phases collapse onto a sin- 
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FIG. 12: (Color online) Parallel self-diffusion coefficient Dy 
as a function of transverse pressure Py . Circles correspond to 
T = 4.0 whereas squares refer to T = 6.0 (see Fig. 3] for color 
code). 

gle curve. In this regime, known as the ballistic regime, 
the MSD increases proportionally to r^. In the ballis- 
tic regime molecules travel freely, that is for sufficiently 
short times interactions with neighboring molecules are 
inconsequential. As a consequence, the specific nature 
of the phase under study does not matter such that the 
MSDs for different phases become indistinguishable and 
can be represented by a unique curve. 

At sufficiently long times, however, the motion of 

molecules becomes diffusive, that is (^Ar^(T)^ oc t. 

In the diffusive regime intermolecular interactions do, 
of course, matter greatly unlike in the ballistic regime. 
Therefore, in the diffusive regime plots in Figs. [TlTa) 
and IllT b) cannot be represented by a single curve but 
differ between I, N, smA, and RN phases. An inspection 
of the MSD for the smA phase reveals that ballistic and 
diffusive regimes are separated by a third, subdiffusivc 

regime in which ^Ar| (r)^ cx r"*-'^-' where 1 < <^{t) ^ 2 

varies continuously. The crossover region of subdiffusivc 
behavior reflects the presence of smectic layers because 

(^Ar| (r)^ is a measure for mass transport in a direction 

normal to the layer plane [sec Fig. [Sjb)]. Diffusive mo- 
tion of molecules out of their original smectic layer into 
a neighboring one is inhibited by the compactness of the 
layered structure. However, within each layer there is suf- 
ficiently little order such that molecules may eventually 
move from their original to a neighboring layer. This 
process demarcates the onset of ordinary diffusive mo- 
tion. Depending on the local compactness of the smectic 
layers it may take some molecules longer than others be- 
fore they exhibit diffusive motion and this is the reason 
why the subdiffusivc regime is characterized by a con- 
stant change in the time dependence of the MSD charac- 
terized by the exponent aij). 

As a quantitative measure of mass transport in the 




Af||(r)/K 



FIG. 13: Root MSD Arx (t) for the average displacement of 
molecules in a direction perpendicular to their long axes as a 
function of the root MSD Ar y (r) for the average displacement 
of molecules in a direction parallel to these axes in units of the 
aspect ratio k. Data correspond to the RN phase at T = 4.0 
and P — 2.6. Solid line is a fit to guide the eye. 

diffusive regime it is convenient to introduce the self- 
diffusion coefficient Z)|| via the expression 

Z?ll= hm ^(Arlir)) (21) 

which corresponds to the long-time slope of the MSDs 
plotted in Figs.[TTJ From that definition and the plots in 
Figs. [TT] one immediately anticipates a self-diffusion con- 
stant which is largest in the RN phase compared with 
the other three phases. In fact, the plot in Fig. W2\ 
shows a dramatic increase of as one enters the RN 
phase whereas lower-pressure (I, N, or smA) phases ex- 
hibit rather small self-diffusivity. In particular, the smA 
phase is characterized by nearly vanishing self-diffusion 
constants which can be rationalized as above where we ar- 
gued that the relatively compact layered structure makes 
it difficult for molecules to diffuse out of their original 
layer and penetrate into a neighboring one. The dra- 
matic increase in mass transport in the direction of n in 
combination with nearly perfect nematic order prompted 
us to refer to liquid crystals in the RN phase as "superne- 
matics" [U. 

More detailed insight into the dynamics of mass trans- 
port in the RN phase can be gained by also considering 
the MSD for molecular displacements in the direction or- 
thogonal to the molecules' long axes [13] to which we refer 
as (Ar^ ('''))( analogy with Eq. ([^0]) where, of course, 

we replace r| by its counterpart r;^ = — {ui ■ ri)ui. 
Plots of the root MSD's for perpendicular versus paral- 
lel molecular displacements Ar^n = ^ (^Ar^ || (r)^ in 

Fig. [13] reveal that dynamically distinct regimes exist. 
In the ballistic regime the plot in Fig. [T2] increases lin- 
early with slope 1 which can be rationalized as follows. 
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Over the relatively short times characteristic of the bal- 
listic regime molecules move freely in space. Hence, their 
self-diffusion is largely controlled by the free volume sur- 
rounding a molecule. In the RN phase the shape of the 
free volume should be closely related to the shape of an 
individual spherocylinder. Therefore, rescaling Ary (r) 
by the aspect ratio k, and expressing Arj^ (r) in units 
of a one anticipates the plot in Fig. [T^] to increase lin- 
early with a slope of 1 in the ballistic regime where both 
Arj^ (r) and Ary (r) are small compared with a. If this 
condition does not hold any longer intermolecular inter- 
actions begin to matter and self-diffusion is no longer 
ballistic. The plot in Fig. [13] shows that during the ini- 
tial stage of non-ballistic self-diffusion the plot of Ar^ (r) 
versus Ar|| (r) levels off and reaches a plateau. During 
this stage molecules diffuse in string-like conformations 
which remain intact, that is the group of molecules form- 
ing a string-like conformation stays together as an entity. 
However, if both Ar_L (r) and Ary (t) are sufhciently 
larger than a the plot in Fig. [T3] increases again linearly 
but with a somewhat smaller slope of about 0.88 com- 
pared with the initial ballistic regime. Because in this 
second linear regime Arj^ (r) is of the order of several 
tr, diffusive motion must involve the exchange of individ- 
ual molecular units between neighboring supramolecular 
"strings" which illustrates the complexity of the dynam- 
ics in the supernematic RN phase. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we elucidate the connection between 
the formation of string-like conformations of elongated 
molecules in the RN liquid-crystalline phase and en- 
hanced self-diffusivity observed earlier by us |2l| in these 
phases. We employ a combination of MC and MD simu- 
lations performed in an isothermal-isobaric and in the mi- 
crocanonical ensemble, respectively. Under conditions of 
confinement a sequence of mcsophases is observed. These 
involve the more conventional I and N phases and, under 
suitable thermodynamic conditions, a smA phase. At 
higher pressures (densities) the more exotic RN phase 
may form which is a nematic phase with surprisingly 
large nematic order and enhanced self-diffusivity in the 
direction of the nematic director n. Quantitatively, we 
express the extent of mass transport in terms of the self- 
diffusion constant D\\. In the RN phase D\\ may exceed 
that characteristic of the I, N, and smA phases by up to 
an order of magnitude. 

Our analysis of the underlying structure reveals that 
the RN phase is characterized by so-called "strings" , that 
is groups of molecules arranged such that their long axes 
and their center-of-mass distance vectors point in the di- 
rection of n with only very little deviation from a perfect 
arrangement. These arrangements can clearly be seen in 
representative "snapshots" of individual configurations. 
On the contrary, string-like conformations are absent as 
far as I. N, or smA phases are concerned. For the smA 



phase this may seem surprising at first on account of the 
formation of distinct layers in the direction of n. How- 
ever, one has to bear in mind that in each individual layer 
the order of center-of-mass arrangements is still rather 
low which apparently prevents the formation of strings 
between adjacent layers. The reader should also note 
that if one considers a single isolated string that remains 

intact indefinitely one would anticipate ^Ary (r)^ cx y/r 

characteristic of single- file diffusion jsO] ■ In single- file dif- 
fusion the time dependence of the MSD is caused by the 
fact that particles cannot pass each other, a situation en- 
countered experimentally for the motion of molecules in 
the narrow spaces existing in materials such as zeolites 
[Slf , single- wall carbon nanotubes Isp , or colloids con- 
fined to one-dimensional channels |53| . However, in the 
present case one needs to realize that we are confronted 
with a totally different physical situation. As the "snap- 
shots" of configurations show many strings exist which 
arc misaligned in the direction in which they diffuse. At 
sufficiently long times individual molecules are exchanged 
between neighboring supramolecular "string-like" units. 
As a result, ordinary Fickian diffusion is the consequence 
in the limit of sufficiently long times. 

The interesting question then becomes: why is this 
diffusive motion enhanced? In fact, based upon plots of 
the cquipotential surface one would assume that even in 
a string-like conformation molecules would "feel" most 
comfortable if individual members in both strings are ar- 
ranged such that neighboring strings are perfectly aligned 
because of the energetically favored side-side arrange- 
ment of pairs of molecules. If this were the case one 
would expect self-diffusion in the direction of the string- 
like conformations to actually be slowed down rather 
than enhanced on account of restoring forces acting in 
the direction of n that pull molecules back to relative ar- 
rangements corresponding to locations of potential min- 
ima. However, this rationale ignores the fact that RN 
phases form at relatively high pressures at which indi- 
vidual strings are "squeezed" together in the direction 
perpendicular to the plane of the solid substrates. That 
strings are actually squeezed in that direction is reflected 
by the local density which shows that in the RN phase one 
additional layer of molecules can be accommodated com- 
pared with either N or smA phases at lower pressures. In 
other words, molecules in neighboring strings are pushed 
towards each other. This happens to an extent such that 
the attractive interactions between any pair of molecules 
pertaining to different strings (and therefore the restor- 
ing forces in the direction of n) are overcome at least 
partly. In other words, the average distance between 
the molecules of a pair decreases with Py such that the 
pair slowly "climbs up" the repulsive part of the fluid- 
fluid interaction potential. Thereby, the net attraction 
between the two molecules is reduced as revealed by the 
shift of the taller maximum oiVlr^) plotted in Fig. [HI As 
a result diffusion in string-like conformations is enhanced 
rather than diminished. 
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This mechanism is similar to one proposed for mass 
transport in zeoUtes. For example, based upon exper- 
iments using quasielastic neutron scattering and corre- 
sponding MD simulations, Borah et al. could demon- 
strate that in NaY zeolites severe confinement causes self- 
diffusivity to become maximum under geometrically fa- 
vorable conditions |22| . This effect known as "levitation" 
is most pronounced if the narrowest void in the zeolite is 
comparable to the optimum distance between a diffusing 
molecule and the pore wall which is, of course, deter- 
mined by the underlying fluid-wall interaction potential. 
If the width of the void is either too small or too large, 
self-diffusivity is lower than under the optimum confine- 
ment conditions. Levitation is quite regularly observed 
in zeolites of various kinds (see Ref. [22 and references 
therein). 



At this point one may ask: Is there any experimen- 
tal evidence supporting the enhanced self-diffusion that 
we observe in our simulations? Though the experimen- 
tal data so far communicated in the literature do not yet 
provide a rigorous confirmation of our theoretical predic- 
tions, some of these data are in reasonable qualitative 
agreement with our simulations. For example, extrapo- 
lating longitudinal relaxation rates reported in Ref. |25| 
from the RN to the N phase yields values markedly be- 
low those in the N phase. If referred to the same T these 
relaxation rates correspond to correlation times that are 
notably shorter in the RN compared with the N phase. 
However, this general interpretation of NMR data re- 
mains speculative as long as translational diffusion has 
not definitely be identified as the process governing the 
observed relaxation. 



Direct evidence for the relation between translational 
diffusion and NMR data can be provided by the pulsed 
field gradient NMR (PFG NMR) technique ^ which 
records molecular displacements typically over a 
range. For example, PFG NMR has been applied to di- 
rectly assess the diffusion tensor upon entering the N 
phase (Hsj . In these studies diffusion in the direction of 
the molecules' long axes was found to increase with in- 
creasing nematic order. By the same technique the diffu- 
sivity of n-alkanes in nanochannels was found to increase 
with increasing orientational order jSGj . These findings 
are in line with our data where enhanced molecular or- 
dering is accompanied by increasing diffusivities in the 
direction of n. The powerful combination of PFG NMR 
with magic angle spinning has recently enabled a notable 
increase in both observation times and gradient pulse 
intensities [s^]. As a consequence PFG NMR diffusion 
measurements became possible beyond the limits of mea- 
sur ability existing so far. This concerns in particular the 
first diffusion measurements with liquid crystals confined 
to nanopores [l^. To stimulate a direct experimental 
verification of our present predictions using these novel 
techniques is the primary purpose of this study. 
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Appendix A: Theoretical background 

From a statistical physical perspective the partition 
function Q of the canonical ensemble is a key quantity as 
far as equilibrium systems are concerned. More specifi- 
cally, it can be demonstrated that with a proper choice 
of phase space variables, Q for linear molecules can be 
written as 



31| 



Q 



1 



A5iV 



X 



N 



(Al) 



where Ath = h/ \/2TTinkBT is the thermal dc Broglic wave 
length. 



Z = 



1 



cxp 



I3U{R,U) dRdU (A2) 



is the configuration integral of a system of N indistin- 
guishable molecules, and /3 = 1/kBT (fee Boltzmann's 
constant, T temperature). The factor 1/2^ in Eq. 
corrects for double counting equivalent configurations. 
These arise because of the head-tail symmetry of the 
molecules, that is the equivalence of orientations speci- 
fied by Ui and —Ut. For the isothermal- isobaric ensemble 
employed in the MC simulations of this work 



A= ^ exp(-/?P||As,o) Q 



(A3) 



(Sx,Sy) 



is the relevant partition function [4^ . It is then straight- 
forward to show that the Gibbsian potential 



(A4) 



is the relevant thermodynamic potential such that ther- 
modynamic equilibrium is attained if Q is minimum at 
any fixed values of A^. Py, T, and s^o- However, it should 
be noted that as far as static properties are concerned 
the specific choice of m or I does not matter (i^ l 

However, as far as the temporal evolution of the sys- 
tem is concerned the specific choice of both m and I 
does matter. To see this let us derive the equations of 
motion that govern the temporal evolution in a system 
of linear molecules. Starting from Hamilton's equation of 
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motion the translation of the molecular centers of mass 
is described by 

r,; = Vi (A5a) 
f 

i), = — , i = l,...,N (A5b) 
m 

where Vi and fi = —VnU are velocity and total force 
acting on the center of mass of molecule i. Similarly, 
one may derive equations of motion for the rotation 
of a molecule starting again from Hamilton's equations 
but this time taking angular velocity oji and Ui as the 
conjugate canonical variables rather than linear momen- 
tum and center-of-mass position as in the derivation of 
Eqs. (|A5|) . The resulting expressions are given by 

Ui = us y. u (A6a) 
^ = i = l,...,iV (A6b) 

where C = Xu) [see Eq. (P)] , 

T, = u, X g,=ux gl (A7) 

is the torque acting on each molecule, and, following Fin- 
cham [5^, Qi = —Vu-U is the so-called "gorque". The 



gorque may be thought of as a turning "force" that effects 
changes in molecular orientation. For linear molecules 
obviously only the gorque perpendicular to the molec- 
ular axis matters. It is obtained through the projection 
of gi onto Ui via 

=9i-{gi-u.t)ui (A8) 

Moreover, in the case of linear molecules the angular 
velocity must always be perpendicular to the molecu- 
lar axis, that is u^i ■ Ui = 0. This makes it convenient 
to replace the angular by the orientational velocity Wi 
such that Eqs. (jA6p can be replaced by the equivalent 
but, from a numerical perspective, much more convenient 
equations 

Ui = Wi (A9a) 
= ^-w^u„ i = l,...,N (A9b) 

The derivation of Eqs. (jA9p requires a little bit of ex- 
tra but straightforward vector algebra [s^. Finally, 
Eqs. (|A5p and (|A9|) show that m and I determine the 
timescale of translational and rotational motion, respec- 
tively. 
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